Scaling and dynamics of sphere and disk impact into granular media 
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Direct measurements of the acceleration of spheres and disks impacting granular media reveal 
simple power law scalings along with complex dynamics which bear the signatures of both fluid and 
solid behavior. The penetration depth scales linearly with impact velocity while the collision dura- 
tion is constant for sufficiently large impact velocity. Both quantities exhibit power law dependence 
on sphere diameter and density, and gravitational acceleration. The acceleration during impact is 
characterized by two jumps: a rapid, velocity dependent increase upon initial contact and a similarly 
sharp, depth dependent decrease as the impacting object comes to rest. Examining the measured 
forces on the sphere in the vicinity of these features leads to a new experimentally based granular 
force model for collision. We discuss our findings in the context of recently proposed phenomenolog- 
ical models that capture qualitative dynamical features of impact but fail both quantitatively and 
in their inability to capture significant acceleration fluctuations that occur during penetration and 
which depend on the impacted material. 
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PACS numbers: 45.70.-n,83.80.Fg,47.50.-d,46.35.+z 

I. INTRODUCTION 

Collisions with complex particulate materials occur in 
diverse situations ranging from asteroid impact [l| to the 
penetration of a running crab's leg into beach sand @. 
Accordingly, collisions with granular media have long 
been investigated [3j|, and like many areas of granular 
research, are being ^actively explored today in experi- 
JTo|. simulation 0, [Tj, and the- 
1131 . [14| . However, because the physics 
of such events must account for both fluid- and solid- 
like behavior during impact, understanding remains lim- 
ited. No comprehensive continuum theory exists for even 
the relatively low impact velocity of a rock dropped into 
beach sand from an outstretched hand. 

Recent experiments and simulations of low velocity 
(below 5 m/sec) impact with granular media have mainly 
considered the depth to which an object penetrates be- 
fore stopping. These studies have investigated how the 
penetration depth scales with various system parameters. 
Durian's group [flll5j performed experiments at low col- 
lision velocities v (maximum penetration depth was ap- 
proximately a sphere diameter) on a variety of spheres 
of differing radii R, sphere density p s , and granular par- 
ticle density p g and found that the penetration depth 

d scaled as d ~ vV 3 { — ) 1 / 2 i? 2 / 3 . de Bruyn's group [j| 



found a different scaling: d ~ ^(^) 1/2 i? 1/2 . The latter 
experiments were conducted at larger impact velocities 
and with spheres of higher density such that maximum 
penetration depths were much greater than the sphere 



diameter d ^> 2R. In two-dimensional (2D) disk simula- 
tions, Tsimring's group found that d ~ vt^ 5 (—) 2 ^ 5 R 3 ^ 5 - 
Phenomenological models have been proposed to ac- 
count for experimental results of penetration depth from 
spheres impacting gran ular materials in a gravitational 
field d, 0, \XL El 111. In all models, the force on a 
vertically falling object impacting a horizontal granular 
medium is written in the general form 



d 2 z 

m— = -mg + F d , 



(1) 



where m is the mass of the impactor, z is the displace- 
ment of the lowest point on the object below the initial 
free surface of the grains, g is the gravitational accelera- 
tion, and Fd is the drag force due to the presence of the 
granular medium. Impact models typically represent the 
drag force as the sum of two terms: 



F d = F z 



(2) 
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where v = dz/dt, and F z is posited to be a 
frictional/hydrostatic term with different proposed 
forms pjj, [l4| depending on the impact regime, and im- 
pactor geometry. For example [l3l |. assumes that an im- 
pacting sphere experiences hydrostatic-like forces, bas- 
ing their arguments on experiments of [16j in which the 
force on a flat intruder moving at a low constant ve- 
locity increases linearly with depth below the free sur- 
face. They argue that F z for shallow impact increases as 
F z = rjp g gz 2 D with r\ a parameter dependent on grain 
properties (e.g. angle of repose), and p g . For deep im- 
pact, the free surface is assumed to move along with the 
sphere such that the bottom of the object is always a 
fixed distance zq below the local free surface; thus F z 
approaches a constant with value F z = rjp g gz 2 D. 
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FIG. 1: Three sequential images (overhead view) of an R = 
1.91 cm steel sphere impacting glass beads at v c = 2.86 m/sec 
with (a) the acceleration and (b) the velocity and penetration 
depth of the sphere during the collision (the three dots (•) 
in (a) correspond with the images). The depth is defined 
as the distance from the lowest point on the sphere to the 
initial free surface of the grains. Also in (a) is a sketch of the 
instrumented projectile (not to scale) showing a single-axis 
accelerometer embedded in a sphere. 



force also includes a term linear in velocity [!, [|| [T2|; 
however, the experiments in 0, Il2| are in a much higher 
velocity regime, 700 m/sec, than the regime examined 
by us and in the other studies cited here where impact 
velocities are typically less than 5 m/sec. 

Since at least 1742 [3| various force laws for granular 
impact have been proposed and their associated pene- 
tration scalings discussed; however, there have been no 
detailed experimental three-dimensional (3D) studies of 
the forces that the colliding object experiences during 
impact. Measuring and understanding the forces exerted 
during impact is clearly important as demonstrated by 
the recent and surprising experimental finding that a disk 
comes to rest in a time t c independent of the initial im- 
pact velocity Fits of position versus time indicated 
that the acceleration during penetration was constant 
with magnitude dependent only upon v c . 

Accordingly, we describe here direct measurements of 
the forces exerted on a sphere during penetration of a 
granular medium. Integrated force measurements show 
how the penetration depth scales in deep impact exper- 
iments (our data most closely follow de Bruyn's scal- 
ing [5|) while a systematic study of collision duration 
reveals simple scaling with system parameters. We ex- 
amine key features of the dynamics using the accelera- 
tion data and then use them to guide an examination of 
the forces at both high and low depths and velocities. 
Constrained by experimental force data we propose an 
equation to describe the forces during granular impact. 
Finally, we show that while our new impact force equa- 
tion and those referenced above account for some of the 
features we observe, a wealth of dynamics associated with 
force fluctuations remains to be understood. 



Durian's group [3| fitted penetration depth scaling 
data from experiment to derive a functional form for the 
frictional drag in the shallow penetration regime such 
that, 

F z = mg + mg[3(z/d ) 2 - 1] exp(-2|z|/d 1 ), (3) 

where do is the penetration depth for an object impact- 
ing the surface with initial collision velocity v c = 0, and 
d\ = m/a. These two parameters are independent of v c . 
Similar to the model of Tsimring and Volfson [l3[, Eq. 
[3] also scales like z 1 for z <§C D, and approaches a con- 
stant (mg) for sufficiently large impact depth. Recently 
Durian's group .17] has proposed that a linear F(z) (sim- 
ilar to the model proposed in [l8[) provides a better fit. 

To summarize, the granular medium is modeled by a 
force law with a hydrodynamic drag term proportional 
to the square of the velocity (which dominates at high 
velocity and thus deep penetration) and a term which 
accounts for a depth dependent static resistance force 
which dominates at low speeds and thus shallow penetra- 
tion depths. Other studies have proposed that the drag 



II. EXPERIMENTAL DETAILS 

Figure [1] shows a sketch of the experimental apparatus. 
An MEMS IC accelerometer (Analog Devices ADXL150) 
with a range of ±50 g's was mounted on an aluminum 
plug which was inserted into a hole or glued to a flat 
on top of the impactor; a small diameter tube isolating 
the accelerometer wires from the granular medium. We 
note that small tube does not affect the impact dynamics 
because the collapsing crater does not contact the tube 
until well after the impactor has stopped (see top panels 
of Fig. [T] [191].) We choose acceleration to be positive in 
the direction opposite gravity. Impactors were dropped 
into the granular medium from heights of 0.01 m to ~ 
2.5 m with corresponding impact velocities v c ranging 
from 0.4 to 7 m/s. 

A variety of impactors and granular materials were 
used, see Tables U and QTJ Since holes were drilled in 
some spheres while flats were made on others and the 
mass of the accelerometer is included in the indicated 
sphere mass, the "effective" sphere density, p s , calcu- 
lated using the masses and radii in the tables does not 
necessarily match the density of the material of which 
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the sphere is composed. Additionally, six 1.3 cm radius 
brass cylinders with different masses were mounted atop 
the 1.9 cm radius nylon sphere, which varied p a from 
1.88 to 9.31 grams/cm 3 . Various containers were also 
employed (see Table llllj) , including a 29 cm diameter by 
40 cm high PVC bucket, a 25 cm diameter by 30 cm 
high aluminum pot, a 50 cm diameter by 75 cm high 
cardboard barrel, and a 10 cm diameter by 10 cm high 
glass jar. Two different procedures were used to prepare 
the granular material: 1) The container was vigorously 
rocked from side to side with decreasing amplitude un- 
til the surface was level. 2) A sieve with outer diame- 
ter approximately equal to the container was placed in 
the bottom of the container, the material poured in, and 
the sieve pulled slowly to the surface. Both procedures 
produced reproducible dynamics. The majority of ex- 
periments were conducted with bronze, steel, or nylon 
spheres or the bronze disk impacting the 0.25-0.42 mm 
glass beads prepared via rocking in the 30 cm diameter 
PVC container filled to a depth of approximately 25 cm. 
Additionally, to vary the effective gravitational acceler- 
ation during impact, an Atwood machine was used to 
drop the bucket containing the granular material with 
accelerations ranging from zero to nearly —g. 



TABLE III: Containers 



Container 


Diameter 


Height 


Wall Thick- 




(cm) 


(cm) 


ness (cm) 


Cardboard Barrel 


50 


75 


0.3 


PVC Bucket 


28.5 


38 


0.3 


Aluminum Pot 


25 


30 


0.33 



III. IMPACT AND DERIVED QUANTITIES 

Figure [1] shows the acceleration, velocity, and posi- 
tion of a sphere colliding with a granular medium. We 
measure the acceleration directly and integrate to ob- 
tain position and velocity. At impact, the acceleration 
increases rapidly as the material suddenly applies force. 
The acceleration decreases as the sphere penetrates into 
the medium and finally comes to rest at a finite pene- 
tration depth d in a time t c . We discuss the details of 
the acceleration profile in Section IV. We begin with a 
discussion of the dependence of d and t c on the initial 
impact velocity v c as well as sphere density p s , sphere 
radius R, and gravitational acceleration g. 



TABLE I: Granular Media 



Material 



Size (mm) pbuik (grams/cm ) r 



Glass spheres 0.25-0.42 1.56 23° 

Aluminum shot lxl 1.62 31° 
(dia.xlen.) 

Millet seed 1.2 0.72 28° 

Bronze spheres 0.05 & 0.17 5.49 24° 



TABLE II: Impactors 



Impactor 


Radius (cm) Mass ( 


grams) 


Steel Sphere 


9.5 


34 


w 


1.3 


83 


w 


1.5 


130 


35 


2.0 


287 


3? 


2.5 


531 


33 


3.5 


1437 


33 


4.0 


2099 


33 


4.5 


3055 


33 


5.0 


4079 


Bronze Sphere 


1.3 


64 


33 


1.9 


201 


35 


2.6 


518 


Nylon Sphere 


1.25 


9 


33 


1.9 


18 


Bronze Disk 


1.0 


10 



(5 mm high) 



A. Penetration depth 

Previous studies 0, [H, [3] have discussed the depen- 
dence of the penetration depth d on the impact velocity 
or, equivalently, the total change in potential energy of 
the impactor (refer to the discussion in the introduction 
for scalings obtained in this earlier work). Figure [2] shows 
that for large enough v c , such that a sphere penetrates 
more than approximately its radius, d increases linearly 
with v c . The linear scaling agrees with the data and the 
scaling proposed in Q. We do not investigate shallow 
impact [14[ in which this scaling is expected to be modi- 
fied [HI due to the varying effective cross section of the 
sphere for z < R. 

We systematically vary R and p s to determine how 
penetration depth changes with these parameters. Since 
our data indicates the dependence of d on v c is close to 
linear for d > R, such that d = do + v c Ad/ Av c , we com- 
pute the slope, Ad/Av c , and intercept, do, of d vs. v c and 
plot these as functions of p s and R as shown in Fig. [3] 
We compare our results for Ad/Av c vs. v c to the scaling 
proposed by de Bruyn d ~ (p s / p g ) 1 l 2 R 1 l 2 v c , and our 
results for d vs. v c to the findings of Ambroso et al. [3, 
do ~ (p s / p g ) 3 / 4 R, since the former reference makes no 
predictions for the scaling of do- As Fig.^a, b) demon- 
strates, the slope of d vs. v c scales with the sphere density 
as d ~ (ps/pg) 1 / 2 and with the sphere radius as d ~ R 1 / 2 
as in Figure Etc, d) indicates that better fits (dashed 
curves) for the scaling of do as a function of p s j p g and R 
are obtained with exponents of 0.59 and 0.75 respectively 
as opposed to the predicted values (solid curves) of 0.75 
and 1. However, comparison of the measured and pre- 
dicted values of the exponents is not strictly valid since 
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Variation of p s 



Variation of R 



v c (m/s) 



FIG. 2: Penetration depth vs. impact velocity for a bronze 
sphere (R = 1.91 cm) impacting glass beads. Solid line is a 
fit with d = (Ad/Av c )v c + do with Ad/Av c = 0.0189 sec and 
do = 0.0185 m. Inset: residuals of the fit. 



our fits apply to the d > R region where d oc v c whereas 
those in [lj] are for d < R and give the actual penetration 
depth of a sphere released at the surface (i.e. v c = 0). 

Combining these scalings, we write the penetration 
depth for d > R as 



iRp, 



d = CiV c * h C 2 {p s /Pg) 

V 9Pg 



0.59^0.75 



(4) 



where C\ and C2 are constants. Figure [4] shows both 
the unsealed penetration depth data and the collapse ob- 
tained using Eq. 2] Although we did not systematically 
investigate the influence of particle (grain) diameter, r, 
on penetration depth, Eq. Q] suggests Ad/Av c ~ r° and 
d ~ rV 4 . 

Figure [5] shows that the penetration depth is affected 
by finite container size. Decreasing the container diame- 
ter by approximately a factor of two decreases the slope 
of d vs. v c by about 1.3 for v c ~ 2 m/s. There is lit- 
tle apparent difference in d for smaller v c . Accordingly, 
the apparent sublinear dependence of d on v c evident in 
the collapsed data in Fig. 0Jc) might result from the fi- 
nite diameter of the container, which, as Fig. [5] indicates, 
becomes more significant the further the impactor pene- 
trates. 



B. Collision duration, t c 

1. Independence of t c on v c 

In a previous quasi-2D study of disks impacting smaller 
disks confined between narrow sidewalls Ciamarra et al. 
found that the collision duration t c was independent of 
v c Q- They attributed this to a constant acceleration 
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FIG. 3: (a, b) Slope and (c, d) intercept of linear fits of 
penetration depth vs. sphere impact velocity in glass beads for 
a 1.9 cm nylon sphere (left column) with p s — 1.88, 2.08, 2.85, 
3.91, 5.03, 7.97, and 9.38 grams/cm 3 (blue, red, black, green, 
cyan, magenta o and blue +) and steel spheres (right column) 
with R = 0.95, 1.27, 1.51, 1.98, 2.46, 3.49, 3.97, 4.52, and 
5.00 cm and corresponding masses m = 34.23, 66.3, 112, 287, 
531, 1437, 2099, 3055 and 4079 grams (blue, red, black, green, 
cyan, magenta o and blue, red, black +). The solid curves in 
(a, b) are proposed scalings from 0] with d — do^ (ps/pg) 1 ^ 2 
and (c, d) from 44) with d ~ (p s /p g ) 3/i R. In (c, d), the 
dashed curves are power law fits with exponents 0.59 and 
0.75 respectively. 



during impact whose magnitude was linearly dependent 
on v c - Our acceleration data for a 3D disk reveals that 
t c (defined here as the time from impact to arrest, see 
e.g. Fig.[T]) is nearly independent of v c for v c ^ 1.5 m/s, 
see Fig. [5] We denote this velocity independent colli- 
sion duration as to. Unlike in [4], but as is the case 
for spheres, we find that the acceleration is not constant 
during the penetration phase (see inset of FigO; we re- 
turn to analysis of the acceleration in Section IV. Below 
v c ~ 1.5 m/sec, t c decreases and as v c — * 0, t c — > 0. 
This decrease at low velocity is a consequence of the fi- 
nite yield stress of granular materials (incorporated in 
the model of de Bruyn [5[), which allows support of a 
distributed finite load without penetration [20( . 

We now make an argument for the dependence of t c on 
v c for disk impact based on the model of [l3j in the deep 
collision regime. Using Eqs. (1) and (2) with constant 
F 7 . appropriate for the regime of deep impact velocity 
[13j . we write 



dt = 



dv 



(F z /m - g) + 



(5) 



We supplement this equation with the boundary con- 
dition for a disk such that v(t = 0) = — v c . Integra- 
tan" 1 a/ (F z —mg)v 
a(F z / m' 2 — g / m ) 



tion yields t = 



C\ with C\ 



5 



Variation of p s Variation of R 




v c (Rg) ~" 2 v c (Rg) 



FIG. 4: Raw data and collapse using the scalings from Fig. [3] 
First column (a,c): Penetration depth and scaled penetration 
depth vs. v c for varying p s . Second column (b,d): Penetra- 
tion depth and scaled penetration depth vs. v c for varying R. 
Symbols are the same as in Fig. [3] 
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1 2 3 4 5 

v c (m/s) 

FIG. 5: Container diameter affects penetration depth. A steel 
sphere with R — 2.5 cm impacting glass beads in cylindrical 
containers with diameters 28 cm (•) and 14.6 cm (o). Both 
containers are filled to a depth of approximately 30 cm. The 
data (•) is reproduced from Fig. [4] The sketches show the 
relative sphere and container sizes. 



tan 1 \J cx.1 (F z —mg)v c 
a{F z /rn 2 ~ g / m) 

We define the collision time t c as when the object 
comes to rest, v = 0. Solving for this time yields 

t = tan' 1 y/ a I (F z ~ mg)v c 
a(F z /m 2 — g/m) ' 

This solution is plotted in Fig. [S] and is a good fit to 
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FIG. 6: Collision duration for impact of an R — 1.0 cm, 
m = 0.01 kg disk with glass beads. Solid line shows fit from 
the model [13| proposed for large impact velocity (deep pen- 
etration) with a = 0.032 kg/m and F z = 0.112 N. At low 
velocity t c increases with v c which is opposite of spheres (see 
Figs. [3 and 0. 



the experimental data; a and F z are fit parameters. As 
v c — > oo, t c approaches t = tt /2J -^p: " l m _ g ) , a constant 

independent of v c . 

The asymptotic form of t c given above allows a predic- 
tion of scaling behavior of to with experimental param- 
eters. Using F z ~ p g gR 3 and a ~ p g R 2 as predicted in 
) to should scale as 

t ^C(^)v^73 (7) 
Pa 

where C is a constant that could depend on impactor 
and grain geometry, density, friction coefficients, normal 
dissipation and other material parameters. 



2. Scaling of to in sphere data 

The experimental result for the disk data showing that 
t c - > to as v c increases is also seen in our sphere data, 
see Fig. [71 Similarly, the acceleration during collision 
is not constant, see inset of Fig. [JJ The large range in 
system parameters varied in the sphere data set allow us 
to examine for the first time how to scales with radius, 
density and gravitational acceleration. 

In contrast to the disk, in which t c — > as v c — > 0, 
t c increases with decreasing impact velocity for spheres; 
compare Figs. [S] and [7J for v c ^ 1.5 m/sec. We attribute 
this difference to the fact that even as v c — > 0, the small 
initial contact area between sphere and grains due to the 
curvature of the sphere always produces local stresses suf- 
ficient for grain bed yielding, which consequently allows 
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FIG. 7: Collision time t c is independent of collision velocity 
v c for sufficiently high v c (R = 2.0 cm steel sphere into glass 
beads). Inset: acceleration vs. time for four impact events 
showing that while t c is independent of v c , the acceleration 
profile is not. 



the sphere to penetrate for a finite time. Thus, a sphere, 
unlike a disk, always penetrates a finite distance into the 
material even with v r = 0; this penetration regime has 
been examined by [EHEj]. 

For sufficiently high v c , the surface of the sphere in 
contact with the granular medium is expected to be es- 
sentially constant for a large fraction of the collision in- 
terval after the initial impact [HI]. We therefore expect 
the proposed scaling of Eq. [7] for disks to be obeyed for 
spheres as well in the high v c regime. As shown in Fig.0 
we find that, to ~ i? 1 / 2 , and to ~ (g otf ) -1 / 2 , where g ett is 
the acceleration of the falling bucket in the Atwood ma- 
chine. For impact at varying sphere density, the figures 
shows that to ~ (p s / Pg) 1 ^ 4 which implies that C in Eq.[7] 
should vary as (p s p ff ) -1 / 4 . 

Combining these three scalings, we obtain 

to~(£i)VV^ (8) 
Pa 

or in terms of the sphere mass to ~ [M/(Rg 2 p g )] 1 ^ 4 . A 
test of this final scaling is shown in Fig. [9] which shows 
that the t c data collapses well for a large range of densi- 
ties and radii. Comparing the experimentally determined 
to scaling to our model shows that Eqs.[5]and[7|both in- 
clude a y/R/g term but have different exponents in the 
density term (1/4 and 1 respectively). 

IV. DYNAMICAL FEATURES OF IMPACT 

While t c and d are characteristic physical quantities 
associated with any impact into deformable material, we 
hypothesize that the scaling that matches experimental 
data can be obtained with a variety of phenomenological 




FIG. 8: Scaling of asymptotic (v c > 1.5 m/sec) penetration 
time to for impact into glass beads as a function of (a) sphere 
density for the R = 1.9 cm nylon sphere, fit shown with to 
in seconds is to = 0.045(p 3 /p 9 ) 1//4 , (b) sphere radius for steel 
spheres, fit shown is to = 0.0507? 1 / 2 , and (c) effective gravita- 
tional acceleration for an R — 1.98 cm steel sphere, fit shown 
is to = 0.075(g rff /(?)- 1/2 . 



models (see references in Section 1). In contrast, examine 
Fig. [10] which shows the acceleration profiles for impact 
of an R = 1.91 cm diameter nylon sphere for five dis- 
tinct impact velocities. As this figure makes clear, while 
current models capture the average physics involved in 
the impact events (we have fit the data to the model of 
[l4T]). they miss much of the detailed physics involved in 
a collision. As v c increases, this model fails to capture 
both significant acceleration fluctuations as well as the 
underlying form of the acceleration. Figure [TTI shows the 
corresponding growth in the relative error between ex- 
perimental data and the model fit as v c increases. 

This model has a second shortcoming: it predicts that 
the scaling lengths do an d d\ in Eq. [3] are independent 
of v c . Least squares fits of a(t) for fixed sphere mass 
M allowing do and d± to vary are shown for different 
v c in Fig. 1101 The dependence of do and d\ on v c is 
shown in Fig. 1121 Although do is roughly independent 



Variation of p s 



Variation of R 




FIG. 9: Test of scaling derived from data in Fig. [8] Collision 
time vs. v c for (a, c) varying sphere density (nylon sphere 
with added masses) and (b, d) varying sphere radius (steel 
spheres). The data in (c, d) are re-scaled data using Eq. p). 
The effective densities of the R = 1.91 cm nylon sphere are 
p s = 1.88, 2.08, 2.85, 3.91, 5.03, 7.97, and 9.38 grams/cm 3 , 
corresponding to blue, red, black, green, cyan, magenta o and 
blue +. The radii of the steel spheres are R = 0.95, 1.27, 
1.51, 1.98, 2.46, 3.49, 3.97, 4.52, and 5.00 cm with associated 
masses ranging from m = 34 to 4079 grams (blue, red, black, 
green, cyan, magenta o and blue, red, black +). 




0.05 0.06 0.07 
Time (sec) 

FIG. 10: (Color online) Acceleration profiles for differ- 
ent impact velocities (v c = 0.47,0.93,1.99,2.88,3.72 m/sec, 
blue, red, green, black, pink) for a 147 gram nylon sphere 
(R=1.9 cm) impacting polydisperse 0.25-0.42 mm glass beads. 
Dynamical features shown and discussed in the text are the 
peak acceleration a pea k, the jump in acceleration as the sphere 
comes to rest, a s to P (see inset), and the time for the object to 
come to rest, t c . The open symbols are experimental data, the 
solid lines are best fits using the model of Ambroso et al fl4| ; 
parameters of the fit are given in the text and in Fig. 1121 The 
model fits are not shown in the inset. 



50 r 



of v c , g?i increases with increasing v c . Since the model 
fails to account for the dependence of do and d± on v c , 
this indicates that additional physics is needed to fully 
characterize impact dynamics. 



Therefore, we now discuss the detailed acceleration 
profile of a sphere as it impacts a granular medium. We 
first describe two robust acceleration features seen for 
a wide range of impactor radii, densities, and material 
types, as well as for a range of granular materials: a peak 
in the acceleration during collision a pcak and a rapid de- 
crease in the acceleration as the object comes to rest, 
Q 3 top, see Fig. 1101 For impact at fixed parameters (i.e. v c , 
R, p g ,p s ) Fig. [l~3l shows that the acceleration profile and 
associated features are reproducible to within approxi- 
mately 5% from run to run. We demonstrate how these 
dynamical features scale with system parameters. Using 
insights gained from the scaling, we deduce an empirical 
force model by examining the experimental data in the 
extremes of high velocity/shallow penetration and low 
velocity/deep (d > R) penetration. 
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FIG. 11: The increasing deviation in average percent differ- 
ence between experimental data and best fits to the model 
of [3 (see Eq. E} from the data set used in Fig. [10] The % 
difference calculation is averaged over the central 70 % of the 
data and model fit. 
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FIG. 12: (a) The scaling lengths do (•) and di (o) in the 
model of [lj] (see Eq. [3]) as determined by best fits to the 
data set from Fig. [10] are not constants as predicted by the 
model but depend on the impact velocity. 




Time (sec) 

FIG. 13: (Color online) The details of the acceleration profiles 
show only small variation (net < 5%) from run to run for six 
separate impacts of a R — 1.91 cm bronze sphere into glass 
with initial velocity 3.65 < v c < 3.67 m/s. 



fluid- like drag [2l|. For fixed v c , a pca k can increase by 
more than 50% if the container is sharply tapped a few 
times before impact. But, as Fig. [T3] indicates, prepara- 
tion by rocking produces consistent results. 




12 3 4 



v (m/s) 
c 

FIG. 14: (Color online) Peak impact acceleration a poa k of a 
sphere as a function of impact velocity into glass beads for (a) 
nylon sphere with effective density 1.88 and 9.38 grams/cm 3 
and R = 1.9 cm, and (b) steel spheres with radii R — 1.3 and 
3.97 cm. Solid lines are fits of a pea k ~ v 2 for v c 2; 1.5 m/sec. 
The insets shows that as v c — » 0, the vf scaling does not hold. 



A. Peak accelerations 

During the collision the acceleration rises to a maxi- 
mum whose magnitude we denote a pca k- At low velocity, 
the peak is not pronounced, but rather is seen as a broad 
maximum. As v c increases, the peak occurs soon after 
the initial contact of the sphere with the grain surface. 
For v c gl 1.5 m/sec the interval between contact and peak 
acceleration is only a few milliseconds. Figure [14] shows 
that in this regime, for all sphere densities and radii, a pea k 
increases approximately like v 2 . a pea k ~ v 2 is in accord 
with all models @, [H, [3| in the high velocity and low 
depth limits, and implies that for these v c impact flu- 
idizes grains sufficiently for the system to display inertial 



At low impact velocity as v c — > 0, a pea k does not ap- 
proach zero, but instead approaches a finite intercept, see 
insets of Fig. [TJ] Figure [TTJ shows that in the low veloc- 
ity limit the peak acceleration occurs at the end of the 
collision just before the sphere comes to rest. We inter- 
pret this as a change from hydrodynamic-like dynamics 
at high velocities to a regime in which resistance to in- 
ertia is no longer the dominant source of drag. This is 
consistent with Fig. [Mja) which shows the lower den- 
sity sphere scaling like v 2 over a larger velocity range 
than higher density sphere; the latter's velocity decreases 
more slowly resulting in deeper penetration where the in- 
fluence of F(z) is no longer negligible. Since high and low 
velocity regimes are dominated by different physics, we 
postpone further discussion of a pca k until the treatment 
of force laws in Section IIV CI 
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B. Stopping acceleration 

For all v c , the sphere does not come to rest gradually, 
but instead suffers an abrupt decrease in acceleration be- 
fore halting, which is reminiscent of a horizontally sliding 
object stopping due to friction. We denote the magnitude 
of the decrease as a atop . This impact feature occurs for 
all spheres, granular materials, and containers we em- 
ployed in our study. As can be seen in the inset of Fig.Q] 
a stop increases with increasing v c . The models of both 
Tsimring and Volfson [l3j], and Ambroso et al. [3] pre- 
dict jumps in the acceleration as the impactor comes to 
rest (see Fig. fTUj) which are attributed to the dominance 
of the depth dependent frictional/hydrostatic drag term, 
F z . While the models predict a(t) instantaneously jumps 
to zero at the end of penetration (Fig. [TO]), we find instead 
that the rate at which a(t) decreases to zero depends on 
v c ; as v c increases, the transition to v — sharpens (see 
inset of Fig. [TO]) . We comment on this in Section HV CI 
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FIG. 15: (Color online) (a) Scaling of a st o P vs. scaled impact 
depth for an R — 1.91 cm nylon sphere with effective densities 
p s = 1.78, 5.03, and 9.39 grams/cm 3 (blue +, red o, black x) 
with corresponding masses m — 51.6, 146, and 272 grams, (b) 
ciatop scaled by sphere radius vs. impact depth for steel spheres 
of radii R = 1.3, 2.5, and 4.0 cm with masses m = 83.3, 531 
and 2099 grams (red, cyan, o and blue, +). 



d is the depth to which the lowest point of the sphere pen- 
etrates the material (see Section IIII Ap . For varying p s 
[Fig. ll5f a)]. we find that a stop is independent of the sphere 
density - at the same depth spheres of different density 
experience the same acceleration (i.e. force oc mass), fur- 
ther evidence for frictional forces dependent primarily on 
geometry dominating the final stages of penetration. a Btop 
increases as the ultimate penetration depth of the sphere 
increases. For d > R a stop increases approximately lin- 
early until d/R w 4.5 . For d < R, a Btop increases more 
rapidly with increasing d (larger slope). We attribute 
the increase in a 3top with d to an increase in the effective 
contact area between sphere and grains as the sphere 
stops with saturation indicating that the lower half of 
the sphere is fully in contact with solidified grains. We 
return briefly to this point in Section IIV CI 

For steel spheres with varying radii a 3top also increases 
with depth and then saturates. The saturation occurs 
at shallower scaled depth and smaller a atop for increasing 
R. Figure [TOTb) shows that before saturation, a Btop is 
independent of R for R varying by more than a factor 
of three and with corresponding masses varying by more 
than an order of magnitude. Although for clarity, the the 
full data set is not shown in Fig. fTST b). the independence 
of a 3top on R and m holds for our entire steel sphere 
data set in which m varies by more than two orders of 
magnitude. Plotting aR/g vs. d collapses the data onto 
a single curve indicating that a 3top varies inversely as the 
sphere radius. 
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FIG. 16: The end of the collision (after a stop ) displays (a) 
overshoot in acceleration which leads to (b) reversal of the col- 
lision velocity and (c) rebound in position for an R = 1.91 cm 
bronze sphere colliding with glass beads at v c — 2.4 m/sec. 
(Temporal axis is zero at start of collision.) 



In Fig. [TO] we examine a atop as a function of penetration 
depth d for varying sphere density p s and radius R, where 



In the final stage of the collision in the vicinity of 
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Qstop, the sphere responds as if it were in contact with 
an elastic-like medium, see Fig. (TBI The sphere does not 
instantaneously come to rest at v = as predicted by 
models, but instead the acceleration decreases rapidly 
(with finite slope) and overshoots (negative acceleration). 
When the rapid decrease in a(t) associated with a stop be- 
gins, the velocity increases through zero, indicating that 
the sphere has reversed direction and is moving upwards; 
for some impact parameters (see curves Fig.fTBTd.e')). the 
velocity and position can oscillate for a few cycles after 
the initial overshoot. Recent work by Durian's group [ItJ 
also observed an oscillation in the velocity of a sphere at 
the end of the collision. They attributed this to displace- 
ment of the bottom of the container, but for our data, as 
seen in Fig. I17f a). the primary overshoot is largely inde- 
pendent of container size and composition for three dif- 
ferent containers: cardboard (blue), PVC (red), and alu- 
minum (green) (see Table [TTTI for container details). Nor 
does the interstitial air diminish the effect as Fig. HHb) 
shows for glass beads. Other evidence pointing to an in- 
trinsic origin is provided by Fig. 117( c) which shows that 
as the sphere density increases the duration of the over- 
shoot decreases and by Fig. 117( d) which shows that with 
the same impactor and container the overshoot varies 
with the granular material. Additionally, the overshoot 
is observed to decrease with increasing sphere radius (see 
Fig.fTTTe)) but to be largely independent of bed depth for 
depths ranging from 8 cm to 20 cm. We speculate that 
when the acceleration reaches a s top, the material sud- 
denly undergoes a solidification transition and the sub- 
sequent dynamics are a result of the sphere oscillating 
within the now elastic-solid as Fig. [TBTc) suggests. This 
picture does not explain why the characteristic overshoot 
time decreases with increasing sphere density. 



Container Effects 



Air Effects 




0.01 0.02 
Time (sec) 

FIG. 17: The dependence of the final stage of the collision (see 
Fig. 1 16 P on system parameters, (a) Collisions averaged over 
interval with 2 < u c < 5 m/s in a cardboard barrel (blue), a 
PVC bucket (red), or an aluminum pot (green) (see text for 
container details), (b) 7? = 2.0 cm sphere with v c = 1.35 m/s 
dropped into a 7 cm deep bed at atmospheric pressure (black) 
and evacuated to less than 50 mTorr (red). Data is an average 
of 6 collisions for each pressure, (c) R = 1.9 cm nylon sphere 
colliding with p s = 1.88, 2.85, 5.03, and 9.38 grams/cm 3 
(blue, red, green, black) Data is an average over all colli- 
sions with 2 < v c < 5 m/s. (d) Collisions with bird seed, 
glass beads, and cut aluminum wire in aluminum pot (blue, 
red, green, black), (e) Steel spheres R — 1.3, 2.0, 3.5, and 
4.5 cm (blue, red, green, black), (f) Layer depths of 8, 10, 
15, and 20 cm (blue, red, green, black). Glass beads are the 
granular medium in (a-c,e,f), an R = 1.9 cm bronze sphere 
is the impactor in (a, d-f), and the plastic bucket holds the 
grains in (c, e, f). 



C. Discussion of force laws 



The preceding results for a pea k and a Btop suggest that 
hydrodynamic forces scaling like v 2 dominate in the 
high velocity/shallow penetration regime, while fric- 
tional/hydrostatic forces are of primary importance as 
v — > 0. To examine these ideas, Figs. [T5] and [T5] show 
the velocity dependence of sphere dynamics at various 
fixed depths during the collision. For shallow depths and 
varying densities, the accelerations are different but the 
forces exerted by the grains, F = m(a — g), collapse onto 
a master curve such that F ~ v 2 (see Fig. [LSa.dh Thus, 
as assumed by jf| 0, EH E2, E3, E2] , a drag force propor- 
tional to the velocity squared, and independent of mass, 
is a good approximation for shallow depths. Addition- 
ally, [E 0, Efl, E3, El, E3| assume the force is inversely 
related to the sphere cross-sectional area. To check this 
scaling we plot in Fig. ITffl a.c) a/g vs. v/y/Rg and F/R 2 
vs. v at shallow depth and for varying radii. Both quan- 
tities fall onto master curves varying like v 2 , indicating 
that indeed F cx R 2 . 

For deeper penetration and decreasing velocity, the 



force and acceleration no longer vary as v 2 , but are in- 
stead linear in velocity with a non-zero offset at v = 0, 
see Figs. QUci) and HStb,d). The slope of F vs. v 
is independent of mass, see Fig HST f). suggesting a low 
Reynolds number fluid-like drag. However, as indicated 
by Fig. ITyT b.d) the linear velocity coefficient of the force 
varies as R 5/2 rather than R as is the case for Newtonian 
fluids. Linear velocity dependence has been proposed by 
[Bl 13, [r3l . although v c was much larger (v c ~ 700 m/sec 
in taE^I) than in our study. For both varying density and 
radii, the extrapolated v = intercept occurs at constant 
acceleration indicating a force dominated by friction. We 
note that at very low velocities (v < 0.25 m/s) a becomes 
constant - this is the regime of a Btop (see Fig. [20T a) for 
example) and is why a stop ^ fi. 

Summarizing the results of Figs. [TBI and 1191 we write 
the drag force Fd exerted by the granular medium on the 
sphere as the empirical force law, 

F d = n{z)mg + C(z)R 5 ^p s ^/gv + a'R 2 v 2 , (9) 
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FIG. 18: (Color online) The acceleration (top panels (a)-(c)) and the force m(a — g) (bottom panels (d)-(f)) of an R — 1.91 cm 
nylon sphere with different effective densities vs. velocity at three distinct scaled depths d/R during impact. Depths are written 
above each column and correspond to (a,d) initial impact, (b,e) penetration of » 1.5 sphere radii, and (c,f) penetration of » 2.5 
sphere radii. Sphere densities p s = 1.88, 2.08, 2.85, 3.91, 5.03, 7.97, and 9.38 grams/cm 3 correspond to blue, red, black, green, 
cyan, magenta o and blue +. 
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FIG. 19: (Color online) The acceleration (top panels (a),(b)) 
and the force m(a — g) (bottom panels (c),(d)) of steel spheres 
of different radii at three scaled depths d/R during impact. 
Depths are written above each column and correspond to (a,b) 
initial impact, and (c,d) penetration of « 2.5 radii. Sphere 
radii are R = 0.95, 1.27, 1.51, 1.98, 2.46, 3.49, 3.97, 4.52 and 
5.00 cm with corresponding masses, m = 34.23, 66.3, 112, 
287, 531, 1437, 2099, 3055 and 4079 grams (blue, red, black, 
green, cyan, magenta o and blue, red, black +). 



where is a depth dependent constant analogous to 
a friction, C(z) is a drag coefficient also dependent on 
depth, and a 1 is a constant independent of depth. We 



consider this equation to be valid in the "steady" collision 
regime before the sidewalls of the crater start to collapse 
onto the impacting sphere (i.e. before the sudden jump 
in a at a atop ). The variation of /j, and C with depth is 
shown in Fig. 1201 The saturation of fJ-(z) with penetration 
depth in the velocity independent term in Eq. [5] is in 
accord with models of [13l. Il4|. However, as we mention 
in our discussion of a stop , this saturation appears to be 
associated with a constant effective contact area between 
the sphere and grains rather than a hydrostatic pressure 
since [i is independent of mass [e.g. see Figs. I18f c) and 
QjOTb)] ■ lu Eq.[9]the linear velocity term follows the model 
in Q for the scaling of Fd with R and g. However, we 
find that the data scale better as p s while || report that 
the scaling goes as ^fpl- 

Finally, despite its agreement with our data in the high 
velocity/shallow depth and low velocity/deep depth, the 
empirical model proposed in Eq. [9] is incomplete: the 
force in the intermediate stages of collision exhibits a 
more complicated non-monotonic dependence on v(t) as 
seen in Fig. |2"T1 



D. Fluctuations and Material Dependence 

For sufficiently low impact velocity, the model of [3| 
captures the shape of the acceleration profile, see Figs. [TO] 
and [TTJ As v c increases, the relative difference be- 
tween the experimental data and the model prediction 
increases. For low velocity, the acceleration is concave 
down. According to the model, the force on the sphere 
increases as the sphere penetrates into the medium, and 
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FIG. 20: (a) The acceleration of an _R = 1.9 cm bronze 
sphere vs. velocity at final stages of collision with glass beads 
(fixed depth d/R = 4.3). The dashed line shows how a atop 
is determined and the solid line shows the linear fit region 
for this fixed depth. Compare to Fig. HSf c - ) and Fig. 118T b). 
(b,c) The coefficients in Eq. [9] vs. penetration depth for im- 
pact of an R = 1.9 cm bronze sphere into glass beads are 
obtained from such fits. The fit region (see Fig. I18[) is for 
0.25 < v < 1.25 m/sec. The fi < region (shaded) in (c) is a 
result of applying a linear fit when the v 2 term is dominant. 




FIG. 21: The acceleration of an R = 1.9 cm bronze sphere vs. 
velocity at intermediate stages of collision with glass beads 
(fixed depth 1.85 < d/R < 1.91) shows non-monotonic be- 
havior (compare to Figs. 1 14l and I15[l . 



exhibit a characteristic structure for higher v c (see also 
Fig. [21]). We attribute the fluctuations in acceleration 
to creation and annihilation of elements of the force net- 
work [22j and are apparently strongly influenced by par- 
ticle shape and also size relative to the impactor. Such 
fluctuations have been observed in many systems but typ- 
ically in a quasi-static regime [23j |. 

As Fig(2U also shows, a pea k increases with the parti- 
cle density p g such that the densest material (aluminum) 
has a peak acceleration approximately six times that of 
the millet seeds (birdseed) at a given v c ; the density ratio 
of the particles is approximately a factor of two (see Ta- 
ble|T|. However, the magnitude of a stop and the overshoot 
at the end of the collision are largest for the millet seed 
and smallest for aluminum which is opposite the behav- 
ior of a pea k- The collision time t c decreases monotonically 
with particle density. 



CONCLUSIONS 



then smoothly levels out to a constant. As v c increases, 
the curves develop upward concavity. Substantial fluctu- 
ations in our experiments appear in the vicinity of the 
change in curvature which implies that they are associ- 
ated with penetration dynamics dominated by the iner- 
tial v 2 term. The models discussed in this paper cannot 
capture such physics as they are purely hydrodynamic. 
In addition, the fluctuations depend on the type of mate- 
rial that the sphere impacts, see Fig. [22] The fluctuations 
are significantly more irregular and occur over shorter 
time scales in the millet seeds (smallest density) and the 
aluminum (largest density) than in the glass beads which 



We have directly measured the forces experienced by 
spheres and disks impacting granular media at collision 
velocities v c < 5 m/sec. By integrating the measured 
acceleration we have deduced scaling relations for pen- 
etration depth and collision time as a function of v c 
and collision impactor parameters and compared them 
to scalings proposed in the literature. We have identi- 
fied robust features of the collision dynamics (a pea k and 
a s to P and described how they scale with v c and intruder 
parameters. These features are in accord with proposed 
models of drag which propose both inertial and frictional 
drag terms. We have shown how the combination of such 
terms determines the surprising finding that time of col- 
lision is independent of velocity. Developing the force 
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FIG. 22: The fluctuations in acceleration during penetration depend on the characteristics of the material being penetrated. 
An R = 1.9 cm, m = 200 gram bronze sphere impacting (a) glass spheres, (b) birdseed, and (c) cut Aluminium wire (see Table 
1 for material properties). Impact velocities (blue, red, green, black magenta) in (a) are v c = 0.45,0.84, 1.82,2.47,3.40 m/sec, 
{b)v c = 0.5, 1.22, 1.74, 2.51, 3.53 m/sec, (c) v c = 0.50, 0.98, 1.79, 2.43, 3.36 m/sec. 



law for penetration empirically from our data, we have 
proposed a new model of the drag force on a sphere dur- 
ing impact which includes constant, linear and quadratic 
terms and have shown how the coefficients in this drag 
relation scale with v c and intruder parameters. We have 
discovered that during impact significant fluctuations in 
acceleration occur which are not described by any exist- 
ing models of impact. Whether these material dependent 
fluctuations fall outside the scope of continuum dynamics 
remains to be determined. 

Our experiments provide another example of the rich 
dynamics and fundamental physics found in the inter- 
action of objects with granular media in the mixed 
fluid/solid regime. While there has been much progress 
in the theory of free granular flow (hydrodynamic 
regime) [H, [2^, [26|, |27[, models describing the mixed 
regime examined here are not rooted in fundamental 
physics but are largely phenomenological. Our data 
can constrain and inform development of models in this 
regime. We hypothesize that to accurately capture the 
dynamics we have observed (including fluctuations and 
stopping acceleration), the physics of rapid fluidization 
and solidification must be included. Beyond discrete sim- 
ulations multi-phase modelling like that proposed in 
[28| might be applicable or perhaps such features will 



have to be captured by statistical models. 

Finally, in addition to probing challenging and un- 
solved problems in the physics of granular media, models 
describing the impact regime we have studied here are 
important to many areas of science and engineering. One 
area with increasing relevance is control of locomotion in 
organisms and robots [29|, l30(. Movement often occurs 
on complex media [3l| and there is a need to understand 
limb interaction with complex substrates [HJ • While the 
scalings of such properties as penetration depth and col- 
lision time might be largely insensitive to intruder and 
grain geometry, the details of the force developed during 
penetration certainly are, and thus theory is needed to 
model such effects. 
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